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Abstract. Until recently our knowledge of the Galactic Bulge stellar populations was based on the study 
of a few low extinction windows. Large photometric and spectroscopic surveys are now underway to map 
large areas of the bulge. They probe several complex structures which are still to be fully characterized as 
well as their links with the inner disc, the thick disc and the inner halo. I will review our current, rapidly 
increasing, knowledge of the bulge stellar populations and the new insight expected towards the Gaia era to 
disentangle the formation history of the Galactic inner regions. 
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1 Introduction 



Two main scenarios, with very different signatures, have been invoked for bulge formation. The first scenario 
corresponds to the gravitational collapse of a primordial gas ( Eggen et al. 1962[ Matteucci fc Brocato| 1990 1 



and/or to the hierarchical merging of subclumps (Noguchi 1999 Aguerri et al. |2001 Bournaud et al. [2009). 
In those cases the bulge formed before the disc and the star formation time-scale was short. The resulting 
stars are old and show enhancements of a elements relative to iron which are characteristic of classical bulges. 
The second scenario is secular evolution of the disc through a bar forming a pseudo-bulge ( Combes et al.|1 990 



Norman et al.fl996 Kormendy fc Kennicutt|[2004| Athanassoula 2005). After the bar formation it heats in the 



vertical direction, giving rise to the typical boxy/peanut aspect. 

Observations of the bulge suffer from the extinction, which varies on very small spatial scales, the crowding, 
and the superposition of different structures along the line of sight. However large scale surveys of the bulge 
are now at reach of the current instrumentation and our knowledge of this complex region is growing very fast. 
We will see here the different characteristics of the bulge as traced by individual stars in our Galaxy and how 
they relate to the different bulge formation scenarios. 



2 Structure 



It is now well established that the Milky Way is a barred Galaxy. Evidence for the presence of a triaxial structure 
in the inner Galaxy came first from gas kinematics, then from infrared luminosity distribution, star counts, red 
clump stars distance indicators and microlensing optical depth. 

But a single structure, with a given semi-major axis and position angle, doesn't seem to reproduce all the 
observations at the same time. The main triaxial structure has a position angle with respect to the Sun-Galactic 
enter direction of about 20° and a semi-major axis around 2.5 kpc. But a thin long bar has been detected within 
10° < I < 27° with an angle of about 45° and a semi-major axis of ~4 kpc (e.g. Cabrera-Lavers et al.|2008 ). The 
main bar dominates within \l\ < 10° with an angle ~ 20° but flattens in the inner regions \l\ < 4° (Nishiyama 



et al.|2005 


Gonzalez et al.|2011a 


(Alard||2001 


Sawada et al.|2004) 



These observations have been found to be well reproduced with a single structure with N-body simulations 
by Martinez- Valpuesta & Gerhard (2011) and Gerhard & Martinez- Valpuesta (2012). In those simulations a 
stellar bar evolved from the disc, and the boxy bulge originated from it through secular evolution and buckling 
instability. They find that the long bar corresponds in fact to the leading ends of the bar in interaction with the 
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Fig. 1. Density profile of the disc particles of the Fux 1999 model projected along |Jf | < 1 kpc (square root density 
scale). The blue straight lines represent the lines of sight along the bulge minor axis b=-4, -6, -8, -12°. It illustrates the 
dominant X-shaped structure at 5° < |6| < 10°. 



adjacent spiral arm heads (see also Romero-Gomez et al.|2011 ). The change in the slope of the model longitude 
profiles in the inner few degrees is caused by a transition from highly elongated to more nearly axisymmetric 
isodensity contours in the inner boxy bulge. We also show that the nuclear star count map derived from this 
simulation snapshot displays a longitudinal asymmetry which could correspond to the suggested secondary 
nuclear bar. 



Looking now at the vertical structure of the bulge, Mc William & Zoccali (2010) and Nataf et al. (2010) 



found two red clump populations coexisting at latitudes |6| > 5° in the 2MASS and OGLE colour-magnitude 
diagrams, corresponding in fact at an X-shaped structure ( Saito et al.|2011 ). This vertical structure is also well 
reproduced in N-body models such as the one of |F^ ( fl999| ) (see Fig. |l[)~or |Li fc S hen] ( [2012] ). 

The main structures of the inner galaxy are therefore now well understood as being shaped by secular 
evolution. As we will see in the next sections, the other characteristics of the inner galaxy are not pointing 
towards this single formation scenario. Those lead Robin et al. (2012) to improve the modelling of the 2MASS 
star counts with the Besancon model by fitting both a triaxial boxy shape (with a slight flare to reproduce the 
double clump) plus a longer and thicker ellipsoid. They note that inner asymmetry is not fitted in this model, 
nor the long thin bar. 



3 Metallicity 



Frogel et al. 


1999 


Zoccali et al. 


2008) 



A metallicity gradient has been detected along the bulge minor axis ( Frogel et al 
In the inner regions (|6| < 4°) no significant gradient in metallicity has been found ( Ramirez et al.||2000 Rich 
etaLl[2007| ). 
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Hill et al. (2011 ) show that the metallicity distribution in Baade's Window (l=l°,b=-4°) can be decomposed 
in two populations of roughly equal sizes: a metal-poor component centred on [Fe/H] = -0.3 dex with a large 
dispersion and a narrow metal-rich component centred on [Fe/H] = +0.3 dex. The same bimodal distribution 
is seen in the sample of 26 microlensed bulge dwarfs of Bensby et al. (2011 ). Their peaks are at [Fe/H] = -0.6 
and +0.3 dex. The discrepancy in the metal-poor peak metallicity should be solved by larger homogeneous 
spectroscopic samples. 

Therefore there seems to be a mix of populations in the bulge, the metal rich population fading while moving 
away from the plane, while those two populations are mixed close to the plane. 

The two red clumps detected in photometry at |6| > 5° has been shown recently to share the same metallicity 
( De Propris et al.|2011 Ness et aT]|2012 ). Only the metal-rich population of the bulge show this split in the red 
clump, implying that the disc from which the boxy bulge grew had relatively few stars with [Fe/H] < -0.5 dex 
( |Ness et al.|[20T2| . 



4 Kinematics 



Babusiaux et al. (2010) show that in Baade's Window the metal rich population presents a vertex deviation 



typical of bar-like kinematics (confirming the results of Soto et al. (2007)), while the metal poor population 
shows isotropic kinematics, close to what can be expected from a classical bulge or a thick disc component. 
They also show that the radial velocity dispersion change with metallicity. The metal poor population shows 
the same velocity dispersion along the bulge minor axis, while the metal rich population goes from bar-like 
high velocity dispersion to disc-like low velocity dispersion while moving away from the galactic plane. This 



corresponds to the kinematic behaviour predicted by the N-body model of Fux ( 1999 ) which contains both a 
disc, which evolved to form a bar, and a spheroid component. However the spheroid component of the Fux 



( 1999 ) model has a too high velocity dispersion compared to the observations. The metal poor component 



observed could correspond to a smaller mass spheroid or a thick disc component. 

The BRAVA (Bulge Radial Velocity Assay) survey ( Kunder et al.|20i2 ), observing at b = —4, —6, —8°, shows 
that the bulge stars follow a cylindrical rotation. Their radial velocity distribution is very well reproduced by 



an N-body model of a pure-disc Galaxy by |Shen et al. (20101, concluding that any classical bulge contribution 
cannot be larger than ~ 8% of the disc mass. However Saha et al. (2012) indicate that the formation of a 



rapidly rotating bar can spin up a pre-existing low-mass classical bulge, leading the classical bulge to develop 
also cylindrical rotation. 



Looking at latitudes b > 12°, Ibata & Gilmore (1995) found that the outer bulge may be kincmatically 



linked to the dissipated core of the halo. 



5 Abundances 



Bulge stars have been shown to be enhanced in a elements (McWilliam & Rich 1994| Zoccali et al. 2006 



Fulbright et al. 2007 Lecureur et al. 2007), suggesting a short formation time-scale. Indeed, in a regime of a 



very fast star formation, most of the stars form with high [a/Fe] ratios, due to the predominant pollution by 
core-collapse supernovae. In slower star formation scenarios, Type la supernovae have the time to contribute, 
decreasing the [a/Fe] ratios. 



However the metal-rich population of the bulge is not enhanced in a elements (Hill ct al. 2011). Gonzalez 
et al.| (2011b) show that those metal rich stars, showing low [a/Fe] ratios at b = -4°, disappear at higher Galactic 
latitudes, in agreement with the observed metallicity gradient in the bulge. They also show that metal-poor stars 
([Fe/H] < -0.2) show a remarkable homogeneity along the bulge minor axis, in agreement with the homogeneity 
observed in their kinematics. 

Several recent studies, based on homogeneous data, also highlight the similarity in term of average metallicity 
and elemental trends between the metal-poor bulge sample and the thick disc ( Melendez et al.|2008 Ryde et al. 
2010[ |Alves-Brito et al.||2010| |Bensby et aT|2011[ ). | Johnson etaL] ( |2012| ) also points out similarities at b = -8° 
between the most metal-poor bulge and halo abundances. 



6 Ages 



The study of bulge colour magnitude diagrams indicate that the bulge is mainly old (>10 Gyr) (Zoccali et al. 
2003 Clarkson et al. 2008). However with infrared photometry van Loon et al. (2003) detected, in addition 
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to the dominant old population, an intermediate age one plus young stars in the inner bulge. |Groenewegen| 
& Blommaertl (120051) found Mira stars with ages 1-3 Gyr up to b=-6°. The study of the density distribution 



of the large-amplitude and long-period variables (Mira variables or OH/IR stars) implies that they trace the 
bar structure of the Galactic bulge (Kouzuma &: Yamaoka 2009). This association of the intermediate age 
population to the boxy bulge component could explain the fact that studies based on those tracers lead to 
larger bar angle (~ 40°, e.g. Sevenster et al. 1999 Groenewegen & Blommaert 2005) than studies based on 



older tracers such as red clump stars (~ 20°, e.g. Stanek et al. 1994 Babusiaux & Gilmore 2005 Nishiyama 



et al. 20051, the old tracers probing the mix of the metal poor component and of the boxy bulge structures. 



Considering the fading of the boxy bulge away from the plane shown previously, one would also expect that this 



intermediate age population represents only a small fraction of the colour-magnitude diagram of Zoccali et al. 

([2003]) at b=-6°. 

However Clarkson et al. ( 2008 ) obtained a proper motion decontaminated colour magnitude diagram with 

, b = -3°) where both the metal poor and the metal rich 
The analysis of the blue straggler stars in this field lead 



a well defined old turn-off in an inner field (1=1 
components would be expected to be well mixed. 



Clarkson et al. (20111 to conclude that the genuinely young (<5 Gyr) population in the bulge must be lower 



than 3%. 



Bensby et al. (2011 1 find that the metal-poor microlensed bulge dwarf stars are predominantly old with ages 
greater than 10 Gyr, while the metal-rich bulge dwarf stars show a wide range of ages, from 3-4 Gyr to 12 
Gyr, with an average around 8 Gyr. The results of Bensby et al. (20111 and the previous discussions about 
the boxy bulge formation seem to be in conflict with the results of Clarkson et al. (2008). To reconcile them, 



Nataf & Gould (2012) suggest an elevated helium enrichment for the bulge relative to that assumed by standard 



isochrones. That would imply that photometric determinations are too old and spectroscopic determinations 
are too young. 



7 Conclusions 



The bulge shows a complex structure both along the major axis and the minor axis. Although the main shape 
of the bulge seems to be driven by secular evolution, the other characteristics of the bulge points toward a mix 
of stellar populations in the bulge. The X-shaped boxy bulge seem to correspond to the metal-rich component 
with solar abundances and bar-like kinematics. A metal-poor component enriched in alpha elements is present, 
corresponding to an older population with a short formation time-scale. This mixed formation scenario has 
been observed in external galaxies (e.g. Prugniel et al. 2001 Peletier et al. 2007 Erwin 2008), as well as 



predicted by several models (e.g. Samland & Gerhard 2003 Nakasato & Nomoto 



2003 Athanassoula 2005 



Rahimi et al. 2010 Tsujimoto & Bekki 2012 1. More observations are currently underway to quantify those 



different populations, their links both in abundances and kinematics, and the links between those populations 
as seen in the bulge and the other populations of the Milky Way: the inner thin disc, the thick disc and the 



inner halo. In the near-infrared, the photometric VISTA VVV survey has made recently its first release (Saito 



et al. 2011). In the optical, Skymapper (Keller et al. 2007) is underway. In spectroscopy, APOGEE in the 



near- infrared ( Allende Prieto et al.|2008 ) and the Gaia-ESO Survey in the optical ( Gilmore et al.|2012 ) are also 
underway. A key piece will of course be added with Gaia, observing 20 million bulge stars (Robin et al. 2005), 
providing parallaxes and proper motion along all the lines of sight towards the bulge and observing red clump 
giants in the bulge in low extinction regions with an accuracy on the paralla x better than 30% , on the proper 
motions better than 1 km s _1 and on the radial velocities around 15 km s _1 (Babusiaux||2011 ). 
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